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Secf/on /. SummarY 

THIS REPORT describes a program for inte- 
grating the restricted three-body problem using 
the TRICE-PB250 system. Difficulties were 
overcome by the use of appropriate variables and 
programming techniques. 

The accuracy of the computer program, together 
with the possibility of continuous visual output, re- 
sulted in a research tool of considerable flexibility 
and versatility. 

Section //. Sfafemenf of the Prob/em 

The restricted three-body problem describes the 
effect of two bodies of finite mass on the dynamic 
behavior of a body (particle) of infinitesimal mass. 
The two bodies of finite mass revolve around one 
another in circular orbits, and the particle of infini- 
tesimal mass moves in their gravitational field. This 
situation is approximately realized, for example, in 
the motion of an artificial satellite in the Earth-Moon 
system. Figure 1 illustrates the geometry of the prob- 
lem (the particle is assumed to move in the Earth- 
Moon plane). 

Section III, Equations of Motion 

In a rotating coordinate system, the motion of the 
particle is described by the following equations (/n = 
the relative mass of the body of smaller mass): 



x + 2y- (1 -ti) 



PARTICLE 



X + M 
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y = y - 2x - (1 -/Lt) 



But major difficulties would be encountered if 
computation were attempted using the equations for 
the rotating coordinate system: 

1. The integration step size would have to be 
changed frequently. 

2. Study of trajectories that pass through a mass 
center would be impossible because these equa- 
tions contain terms with r* in the denominator. 




Figure I . Geometric representation of the restricted 
three-body problem in a rotating coordinate 
system. 



3. The position and velocity terms would vary 
greatly in range. 

4. Scaling would be extremely difficult. 
Therefore, at the suggestion of Dr. Arenstorf, staff 

mathematician of the Computation Division, a Thiele 
transformation was used to overcome many of the 
computational difficulties and also to provide a more 
flexible form of the equations. Use of the Thiele 
transformation provided a further advantage, in that 
broader studies of the problem were possible. 

The Thiele transformation is accomplished by 
letting 



X - /g - M + Vj cos u cosh v 
y = - /, sin u sinh v 



The new independent variable is then defined by 
t 
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FIG. 2. TRICE PROGRAM for the restricted 



Making these substitutions into the rotating form 
results in the following:* 



u" = Vj (cosh* V - cos* u) v' 
+ sia u 

- — cosh V (cosh* V - 3 cos* u) 

lb 



'°"x A^<7*\ 



+ ■— cos u cos 2u 
16 



V" = 



- r >, ( cosh* V - cos* U) U' 



+ sinh V 



* V2 16/ 



cosh V 



j^ cos u (cos* u - 3 cosh* v) 



+ ^ cosh V cosh 2v 



*A prime (') after a quantity denotes the first derivative 
of the quantity with respect to s; a double prime (") de- 
notes the second derivative with respect to s. 
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three-body problem in a rotating coordinate system. 



where ^ = 1 — 2/j,, and h (Jacobi's integral) is de- 
fined by the following equation: 

(u')2+ (v')2 -Ha COSU -coshv 

- - (cosh V - cos u) - — sin^ p (eos u + cosh v)^ 

- rr sinh^ v (o- cos u + cosh v) ^ = 
16 

Via a third transformation, the motion of the 
particle was also described in a star-fixed coordinate 
system. The equation 



t = >4 y (cosh^ V - cos^ u) ds 



was used to generate a new independent variable, 
and the following new equations were introduced: 

Zj = X cos t - y sin t ^2 ~ ^^ ^°^ t + x sin t 



Section IV. The Computer Program 

The configuration of TRICE elements that mech- 
anize the Thiele form of the equations of motion is 
shown on Figure 2. The initial conditions for each 



trajectory were entered into the PB250 in terms of 
the rotating coordinate system variables /*, x^, y^, 

Xo and yo 

The PB250, under control of an interpretive pro- 
gram language, computed a set of initial conditions 
that satisfied the Thiele form of the equations. 

The following calculations were made: 



2. r2=Y(X(, + M-l)2 + y,2 
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V, (cosh^ Vjj - cos^ U|j) 



3. cosh Vg = Tj + Fj 



4, cos u. = rj - x^ 



5. (cosh Vo)^ = (ri + rj)^ 

6. (cosuo)^= (ri - rs)^ 

7. To = V* (cosh^ Vo - cos* po) 

8. sinh Vo = '^cosh Vq - 1' 



9, sin 



uo='/r 



cos Un 



10. CT = 1 - 2;i 



11. 0^2= (i _ 2;^)' 



12, uq' - /j (-yj COS Uo sinh vo + Xo sin Uq cosh Vj) 

13. Vo' = /j ( -yo sin Uo cosh Vo + Xo cos Uj sinh vj) 



15. cr cosh V. 



16. <J cosh V + cos u„ 




17. (7 cos u + cosh Vp 



18. a^ cos u +cr cosh v. 



19. 2 cos* u - 1 





The results of these calculations were stored in 
the registers of the TRICE modules. 

The major function of the TRICE program was to 
integrate the Thiele equations. Analog plots of the 
solution were made in the rotating coordinate sys- 
tem and in the star-fixed coordinate system. 

The Thiele variables were transformed into the 
rotating coordinate system by TRICE modules gen- 
erating the following functions: 



x = %-/^-^V, 
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FIG. 3. TRICE PROGRAM for the star-fixed 
coordinate system. 



The TRICE program for the star fixed transfor- 
mation is shown in Figure 3. 

Section V. Cone/usions 

A. ADVANTAGES OF THE THIELE FORM 

Certain advantages of the Thiele form of the equa- 
tions of motion are worth noting: 

1. The choice of the independent variable 

t 



'V J r^TiJ 



V 

eliminated the need to vary the integration step 
size as the trajectory passed through different re- 
gions of space. An equivalent effect was accomplished 
continuously during the solution as a result of the 
mathematical nature of the Thiele form. This effect 
can be clearly observed when trajectories approach 
either mass. Motion of the particle becomes very 
slow during the approach and increases in speed 
during departure from the region of a mass center. 

2. As was noted earlier, because of terms containing 
1/r^, the study of trajectories that pass through a mass 
center would be impossible using the equations for 
the rotating coordinate system. These terms do not 
appear in the Thiele form. 

3. The terms in the Thiele equations are better 
behaved than the position and velocity terms in the 
rotating set. The sinusoidal terms are bounded be- 
tween -|-1 and — 1 for all real arguments. The hyper- 
bolic functions, in contrast to the position and velocity 
terms in the rotating set, do not vary greatly in range. 

4. The behavior of the new Thiele variables elimi- 
nated many scaling problems. 

B. ADVANTAGES OF THE TRICE-PB250 SYSTEM 

Use of the TRICE-PB250 system, which combines 
advantages of both digital and analog computers, 
contributed to an effective and efficient study of the 
three-body problem. Some of the advantages of the 
system include: 

1. Parallel integration without the problem of drift 

2. Generation of hyperbolic functions without stability 
problems 

3. Continuous visual output provided by the converter 
and plotter 

4. High accuracy and repeatability (6 decimal digits) 

5. Automation provided by the PB250. 

The orbits shown on Figures 4 through 6 were ob- 
tained using the TRICE-PB250 system. The study re- 
quired that the initial conditions be adjusted until 
the trajectory became periodic. Then, using the visual 
output provided by the plotter, a human operator 
determined new initial conditions on the basis of 
symmetry and orthogonality (a periodic orbit is 
symmetrical about the horizontal axis and makes two 
orthogonal crossings). An automatic iteration process 
for finding proper initial conditions appeared too com- 
plex for expression in a useable mathematical form. 
Because of the relatively low cost and the ready 
availability of the system, the operator was able to 




FIG. 4. TWO-LOOP periodic orbit computed using 
the TRICE-PB250, with initial conditions: h = — 

1.197; xorr: 0.994; yo = — 2.ll4;;a = 0.1227747; 




FIG. 5. THREE-LOOP periodic orbit computed 
using the TRICE-PB250, with initial conditions: 
h = — 1.367; xo = 0.994; yo = — 2.0318; fi. = 
0.12277471; m/k = 2/3 " 




FIG. 6. FIVE-LOOP periodic orbit computed using 
the TRICE-PB250, with initial condition: h = — 

45076; xo = 0.9927; yo = — 2.64215; ^ = 
0.122775; m/k = 3/7 retrograde. 



FIG. 7. SEVEN-LOOP 

periodic orbit using 
TRICE-PB250 with in- 
itial conditions: Xo = 

I.OI; Xo = 0; yo =: 0; 

yo = —1.369; 
/82. 




FIG. 8. FIVE-LOOP periodic orbit (per- 
turbed ellipses) in a star-fixed coordinate 
system. 



make the detailed studies required for this type of orbit. 

The orbit shown on Figure 7 is a solution plotted in the rotating 
coordinate system. Figure 8 shows the trajectory as it appears in the 
star-fixed coordinate system. 

Running times for a periodic solution depended on the initial condi- 
tions. Running times as long as 30 minutes were required for very-low- 
energy orbits, and as short as one minute for very-high-energy orbits. 

Results to date have been of significant value as guidelines for further 
theoretical studies. 
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